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The Sudbury Neutrino Observatory (SNO) has given precision data on the flux of the 
®B solar neutrinos and on its flavor content, and these data have provided solid evidence 
for neutrino oscillations PlEl- The SNO experiments measure the yield from the charged 
current (CC) reaction, Uf. + d ^ e~ + p + p, which occurs only for the electron neutrino, 
and the yield from the neutral current (NC) reaction, Ux + d^Ux+p + n (x = e, /x, r), 
which takes place with the same cross section for any neutrino flavor x. In interpreting 
the SNO data the theoretical estimates of the u-d cross sections play an important role. 
Detailed studies of the u-d reaction in the solar neutrino energy region have been carried 
out using the standard nuclear physics approach (SNPA) H], pion-less effective field 
theory based on the power divergence subtraction scheme |S1 E] , and EFT* (or MEEFT) 
[Zj, which is a "pragmatic" version of chiral perturbation theory in the Weinberg counting 
scheme |H] . Apart from the radiative corrections, the uncertainties in the calculated values 
of the u-d cross sections are estimated to be 1%, see Ref. 

The above-mentioned calculations use the weak interaction Hamiltonian dictated by 
the standard model, in which the neutrinos interact with the deuteron only by exchanging 
the W and Z bosons. Meanwhile, if the neutrino has a magnetic moment (let it be denoted 
by iJ,i,), it can interact with the deuteron via a photon exchange as well. This additional 
interaction, contributing to the u^d — >■ UxPn reaction, can affect the interpretation of the 
NC data of SNO at some level of precision. In the following, the electromagnetic neutrino- 
deuteron interaction due to shall be simply referred to as the /i,^-interaction. Since the 
contribution of the /i,y-interaction to the u^d — ^ i^xPn reaction and the NC contribution 
do not interfere with each other (see below), we can decompose the cross section as 



where ctem stands for the contribution of the /i,y-interaction. Akhmedov and Berezin 
(AB) P estimated aEM, using the effective range approximation and the impulse approx- 
imation. Recently, Grifols, Masso and Mohanty (GMM) lOj gave a new estimate of the 
/ij^-interaction effect based on the equivalent photon approximation. GMM reported a 
significant difference between their value of aEM and the value obtained by AB P, but 
the origin of this discrepancy was left undiscussed. 

We carry out here a calculation of aEM which is free from the effective-range approx- 
imation, the impulse approximation and the equivalent photon approximation. To this 
end, we make an appropriate extension of SNPA used in Refs.[SlllI and calculate aEM, 
including both the impulse and exchange-current terms, and with the use of the two- 
nucleon wave functions generated from a high-precision realistic NN potential. We shall 
also briefly discuss some consequences of our calculation in the context of placing upper 
limits to the neutrino magnetic moments. 

We are concerned here with the reaction 



at low energies, where a and b stand for neutrino flavors. The relevant transition amplitude 
tfi consists of two terms, one arising from the standard NC interaction and the other from 



a{uxd u^pn) = aEM + (^nc, 



(1) 




(2) 
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the yUy-interaction. Thus 

, Gp -NCiNCtJ. I 1 ■EMiEM,fi /Q\ 

tfi = ' ' 

where qii = k^ — k'^. The lepton neutral current matrix element, Z^'^''^, is diagonal in the 
weak-flavor eigen-states and is given as 

F = Ub{k')Y{l-J,)Ua{k)Sb,a, (4) 

whereas the neutrino electromagnetic current is given in terms of the magnetic dipole 
moment fiba and the electric dipole moment eta, 

P^./^ = ui,{k')ta^''q,if^,a-etal,)ua{k) , (5) 

where b and a stand for the mass eigen-states of the neutrinos. In Eq.©, j^*" and 
represent the nuclear matrix matrix elements: 

^ = <p{p[)n{p',)\J^\d{P)> {X = NC,EM), (6) 

where J^*" is the weak NC current, and J^*^ is the electromagnetic current. The explicit 
forms of J^^ and J^^^ were described in detail in Refs. [31 13]. At the solar neutrino 
energies it is safe to apply non-relativistic approximation to the one-nucleon (impulse 
approximation) current. The two-nucleon exchange current is derived using the one- 
meson-exchange SNPA model [SI E] . We remark that this SNPA model describes very 
well the n + p ^ d + 'j reaction at low energies, and that, for En < 0.1 MeV, the magnetic 
dipole (Ml) transition dominates, whereas for higher energies the electric dipole (El) 
transition dominates; see Ref. jTTj for details. 

If the neutrino mass is neglected in Eqs.(jH) and (jSl), then l^'^^t^ conserves the neutrino 
helicity while /^^^-/^ flips the helicity. This means that to a very good approximation there 
is no interference between the NC and EM contributions, and therefore the total cross 
section for the u^^d — UxPn reaction can be written as the incoherent sum of the NC 
and EM terms as in Eq.(^. An explicit expression for aNc and its elaborate numerical 
calculation can be found in Ref.pj; our goal here is to evaluate aEAi to a comparable 
degree of accuracy. It is convenient to express a em as an integral over the final two- 
nucleon relative kinetic energy, T^^ = p'"^ /M, where p' = {p'l — v'^tjl, 

rTmax dCEM 
C EM = / dT]y]y — — . (7) 
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where q is the momentum transferred to the two-nucleon system; the neutrino scattering 
angle 6 is defined by k k' = kk' cos 9 in the u-d center-of-mass system, fi^a is the effective 
neutrino magnetic moment for the incoming neutrino of flavor a fT^ IT^ 

if^ua)^ = '^\fiba + (^ba\^ ■ (9) 
b 

In Eq.®, <T^> represents the nuclear reduced matrix element of a multipole operator 
Tx of rank J: 

<Ti> =<{Lf,Sf)Jf\\Ti\\d;J, = l> , (10) 

where X = C, M and E correspond to the Coulomb, magnetic and electric multipole 
operators, respectively. The explicit expressions for < TjJ > can be found in Eqs. (57) 
and (38)-(40) in Ref. |3j. In calculating aEM at the solar neutrino energies, we need to pay 
particular attention to the sizeable flnal two-nucleon interactions. Here, as in Ref.|3], we 
use the flnal two-nucleon distorted S- and P-waves obtained with the use of the Argonne 
V18 potential dH. 

In order to illustrate the possible /ij^-interaction effects on a{vxd i^xpn), we need 
to specify a "typical" value of fi,ya- The radiative corrections within the standard model 
leads to an estimate /i,^ ~ 3 x 10~^^m^/leV[fiB], where fiB is the Bohr magneton [T^IT^ : 
the value of fi^, larger than this estimate would be an indication of new physics. The 
current upper limits obtained from reactor neutrino experiments are /ip^ < 1.0 x 10~^^ fiB 
[TH] and /ip^ < 1.3 x lO^^^/i^ [ini- The limit found using the and from 7r+ and /i"*" 
decays is /i^ < 6.8 x 10~^°/zb jUj. Meanwhile, the limit for the tau neutrino is reported 
to be /It < 3.9 X 10^^/iB[TH].^ As for the astrophysical constraints, it is known that the 
avoidance of excessive stellar cooling requires /i^a < (0.3 — 1) x 10~^^/is (see, e.g., Ref. [T^ ) . 
For the sake of deflniteness, we use here fi^a = 10~^°/i_B as a representative value. This 
completes the speciflcation of all the quantities appearing in Eqs.©, 

Our numerical results for aEM are shown in table 1 and Fig. [T] In Fig. [H the dotted 
curve represents the contribution of the Coulomb and El transitions leading to the ^Pj 
flnal two-nucleon state, while the dash-dotted curve gives the contribution of the Ml 
transition leading to the ^5*0 state. The solid curve shows the value of ctem obtained 
by adding the contributions of all the multipoles in Eq.®. We note that, except for 
the lowest energy region {E^ < 2.8 MeV), aEM is dominated by the Coulomb and El 
transitions to the flnal P-wave states. This should be contrasted with the case of the 
NC reaction, where the Gamow- Teller transition to the flnal ^5*0 state gives a dominant 
amplitude. 

Fig. 121 shows the ratio R = aEM/c^NC calculated from aEM obtained in the present 
work and (Tatc* given in Ref. 1^. R increases with the neutrino energy E^ and reaches 

"'^When neutrino oscillations are taken into account, a limit deduced from scattering experiments should 
be understood as a limit on /i/ = | Ufkfijke~'^^''^\'^ , where / is the flavor of the beam neutrino, k 
and j denote mass eigenstates, and U is the mixing matrix [T^. 
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~ 5 X 10^^ at Ejy = 10 MeV, which is however still very small. The smallness of R 
reflects the fact that in the u-d disintegration reaction, the NC contribution comes from 
an allowed transition, while the EM contribution arises from a 'first-forbidden' transition. 
The latter is intrinsically suppressed by a small factor q <r >, where < r > is a typical 
nuclear size. We remark that the contribution of the /ijy-interaction is much larger for 
elastic u-d scattering than for the u-d breakup reaction, since the elastic scattering can 
occur via an allowed-type transition; it is unfortunate that u-d elastic scattering is at 
present practically impossible to detect. 

To facilitate comparison of our present results with those in the literature, we let 
c"£;a/(TNSKM), (T£;m(AB), aEMiGMM) denote the values of aEM obtained by us, by AB 
ini and by GMM (TU], respectively. To compare (T£;m(TNSKM) and cr£;j\/(AB), we read off 
the numerical value of (T£;m(AB) from Fig.l in Ref. [H] and normalize it to the case of /Xj, = 
W~^^Hb, the value used in our present calculation. Comparison after this normalization 
indicates that cr£;Af (TNSKM) is about 4 times larger than (Tem(AB). As mentioned, our 
present calculation based on SNPA is free from the effective-range approximation and 
the impulse approximation, whereas the calculation in Ref. [D] does involve these two 
approximations. However, this difference in formalisms per se is not expected to lead to a 
discrepancy as large as a factor of ~ 4. For further comparison, we find it illuminating to 
simplify our full calculation by introducing the effective-range approximation and impulse 
approximation; the resulting simplified formula for the El contribution is found to be 
larger than the El part of Eq.(18) in Ref. ^ by a factor of 2. It is likely that the remaining 
discrepancy by a factor of ~2 is of the numerical origin. We have done a little numerical 
exercise of calculating "our version" of (T£;m(AB) directly from Eqs.(18) and (Al) of AB. 
Let ctemIAB; direct) stand for the result of this exercise; the value read off from Fig.l in 
AB is denoted by (T£;a/(AB; figure). If we use ct£;m(AB; direct) instead of cr£;A/(AB; figure), 
then the difference between our result and that of AB becomes a factor of ~2 (instead of 
~4), which is consistent with the above-mentioned difference in the analytic expressions. 

As far as comparison with GMM |10j is concerned, we first discuss <aEM> defined as 



where a = EM,NC and f{Ejj) is the normalized ^B neutrino spectrum. If we use the 
value given in GMM, we find <ai5Af(TNSKM)> ^ <aEMiGMM)>, but that the former is 
smaller than the latter by about 10 %. To gain some insight into this difference, we have 
examined the analytic expression for the cross section in GMM and have confirmed that 
Eq.(16) in GMM is indeed valid in the long wavelength approximation and in a regime 
where the dominance of the electric and Coulomb dipole transitions holds. Furthermore, 
if we calculate aEM directly from Eqs.(16) and (17) of GMM, the result agrees with 
cbm(TNSKM) within 1%. Therefore, the above-mentioned ~10 % discrepancy does not 
seem to be due to the equivalent-photon approximation used by GMM. We have not been 
able to identify the origin of the ~10 % difference found in the comparison of <(Jem>-'^ 

^In GMM's article ^Uli there is the statement that (jEKiiGMM) is smaller than (7£;m(AB) by a 




(11) 
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In these circumstances it is probably warranted to repeat that, if one apphes the same 
method (SNPA) as used here to the calculation of the np 'yd reaction cross section 
[llj . the results agree very well with the experimental data 20 for both the Ml and 
El contributions. Furthermore, the calculation of the cross sections, a{Vf.d — > e~pp) and 
(j{vxd — >■ UxPn), with the use of SNPA has been done by several groups (with different 
degrees of sophistication), and the results are known to be consistent among themselves 
at the level of precision relevant to the discussion of ctem here j2I]. Thus it is our behef 
that the SNPA method used in the present work has been well tested in both formal and 
numerical aspects. 

We now briefly discuss consequences of our calculation in the context of obtaining 
information about the neutrino magnetic moments. In considering the SNO data, we 
find it useful to follow the argument of GMM JU]- The total "NC" deuteron dissociation 
events observed at SNO is the sum of the standard NC events and those due to the 
/x,y-interaction: 

N'^P = NdT<^suN [<aNC> + <(Tem>] , (12) 

where A^^ is the number of target deuterons and T is the exposure time; ^sun is the 
total solar neutrino flux. The empirical *B neutrino flux, ^sno, obtained at SNO 
assumes that the total dissociation events arise solely from the standard NC interaction. 
Thus (^sNO = N'^'^P/lNdT <aNc>], which implies (^sno = ^suNi^ + with 5 =< 
<^EM> I <<^NC>- To extract information on /i^, GMM employed oem = c^ba/IGMM) 
and the value of a^^c obtained by Nakamura et al.|l], and furthermore they assumed 
^SVN = ^SSM, where ^ssm is the *B solar neutrino flux predicted by the standard 
solar model GMM then arrived at ^sno = ^ssm{1 + 6.06 x 10^^/2^), where p, 

is the operational neutrino magnetic moment defined by /t^ = /igi + + fAs- The 
use of ^sNO = (4.90 ± 0.37) x lO^cm-^s'^, and ^ssm = (5.87 ± 0.44) x lO^cm-^s"! 
leads to /i^ = (—2.76 ± 1.46) x 10~^^, where the errors in ^sno and ^ssm have been 
added quadratically. Inflating the error up to 1.96cr, GMM deduced the upper limit 
< 3.71 X lO^^fiB (95% CL). If we repeat a similar analysis using (T£;Af(TNSKM), we 
arrive at ^sno = ^ssm{^ + 5.46 x lO^^/i^). We then should obtain practically the same 
conclusion on as in GMM. This implies that, if the upper limit ~ 10~^°/iB to fii, is 
adopted, the interpretation of the existing SNO data is not likely to be influenced by the 
possible existence of fi,^. 

Recently, an extensive global analysis of the totality of the data from the Super- 
Kamiokande, SNO, and KamLAND has been carried out, with the possible existence 

factor of ~ 2. A comment is in order on this point. As mentioned, (T£;m(AB) in fact can represent either 
o'_EM(AB;figure) or cr £;M(AB;direct). In the former case, we are led to cr_EM(TNSKM) « 4f7_EM(AB;figure) 
« 8aEM{GMM). In the latter case, we are led to cr£;A/(TNSKM) « 2cr£;M(AB;direct) « 4cri5Af(GMM). 
However, as mentioned in the text, we find (T£;m(TNSKM) « aEM{GMM). This puzzling situation can 
be avoided if the statement (tem(GMM)« 0.5(tem(AB) in Ref.pT)] is changed into the statement (in 
our notation) (Tbm(GMM)r:! 2(T£;M(AB;direct). In a private communication Dr. S. Mohanty has kindly 
informed us that he agrees with this change. 
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of jj^i, taken into account [22j. This analysis reports an upper limit /Xg// < 1.1 x 10^^ 
(90 % CL). The value of aEM described in our present communication has been used as 
input in the global analysis in Ref. |23j. 

It is known that a non-zero value of fii, can affect the shape of the recoil electron 
spectrum in electron-neutrino scattering. A similar effect is expected in the u-d elastic 
scattering (which is however at present practically impossible to measure). For the p-d 
disintegration reaction, the effect of the interaction is much smaller than for the elastic 
scattering, but it might be of some (academic) interest to discuss to what extent one can 
distinguish the NC and /ij^-interaction contributions to the neutron energy spectrum in the 
u-d breakup reaction. As an example, we consider the T^n dependences of daEAi/dT^N 
and daNc/dTNN, where T/vat is the final two-nucleon relative kinetic energy. Figl^l gives 
these differential cross sections calculated for E^, = 10 MeV, and for /Xe// = 10~^/xb,3 x 
10~*/iB and 10~^ fiB- These values of fieff are much larger than the upper limits known 
for z/g and z/^, but within the upper limit pertaining to u-r- In FiglHl the NC contribution 
is shown by the solid line, while the /i^ contribution is given by the dotted, dashed and 
dot-dashed curves corresponding to fiejf = 10~^yUB,3 x lO^^/i^ and 10~^ fiB, respectively. 
The NC contribution has a sharp peak in the low T^vat energy region, owing to the final 
^5*0 state dominance for the Gamow- Teller transition. By contrast, the /ij, contribution 
exhibits a maximum away from threshold, reflecting the fact that /ij^-interaction leads to 
the final P-wave states. As illustrated, fieff of the order of 10~'^fiB can affect the shape 
of the energy spectrum of the recoil nucleon in the u-d breakup reaction. 

To summarize, we have studied the effect of the neutrino magnetic moment on the 
neutrino-deuteron breakup reaction, using a method called the standard nuclear physics 
approach (SNPA), which has been well tested for a number of electroweak processes 
involving the deuteron. The present calculation is free from the various approximations 
(mentioned earlier in the text) that were used in the previous estimations. 
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Figure 1: The cross section, aEMi for the v^d i^xPn reaction due to the /Xi^-inter action. 
(Tem has been calculated for /^e// = 10~^°A*s, and is given in units of cm^. 
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Figure 2: The ratio R = ctem/o'nc calculated for /Xe// = 10 as a function of the 

incoming neutrino energy. 
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Figure 3: The neutrino dcuteron disintegration cross sections, a^c and a em, as functions 
of TjvAT- See the text for details. 
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Table 1: The cross section, oem, for the Uxd i^xPn reaction due to the interaction. 
aEM has been calculated for /Xg// = IO^^^/^b, and is given in units of cm^. 
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5.31 


(-49) 


9.2 


4.67 


(-48) 


12.8 


1.45 


(-47) 


16.4 


3.06 


(-47) 


20.0 


5.31 


(-47) 


5.8 


6.36 


(-49) 


9.4 


5.06 


(-48) 


13.0 


1.52 


(-47) 


16.6 


3.17 


(-47) 









11 



